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TER INTRODUCTION 


The simulation of flight conditions using ground based 
facilities has been an essential part of all aircraft 
development. Figure 1 illustrates the various types of 
facilities which have evolved to simulate flight at 
progressively higher speeds and altitudes. Increasing 
energy requirements led to blow-down wind tunnels. The need 
to avoid condensation in nozzle expansions required heaters. 
The very high stagnation enthalpies associated with 
hypersonic flight led to very short duration, intermittent 
techniques. Each of the types shown in Figure 1 was іп 
operation more than 20 years ago. А renewed interest in 
hypersonic simulation has been generated by the National 
Aerospace Plane(NASP) program. 

A recent review of facility requirements for hypersonic 
propulsion system testing by Dunn and Lori [Ref. 1] 
concludes that a single ground based facility will not be 
capable of providing the necessary testing to cover the NASP 
flight envelope. The planned use of Computational Fluid 
Dynamic (CFD) codes to bridge the gap left by inadequate 
ground based facilities will make CFD codes an essential 
part of the design process but, before confidence can be 
placed in these codes, fundamental experiments are needed to 


validate the models incorporated for real gas behavior. 
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Figure 1. Ground Based Testing Facilities [Ref. 1] 
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Recognizing the need for CFD validation, the feasibility 
of establishing a shock tunnel capability by modifying the 
existing Naval Postgraduate School (NPS) three inch diameter 
shock tube, was examined in the present study. The time 
duration and magnitude of constant pressure available in the 
reflected region of the shock tube were evaluated. The 
concept of tailoring was examined using two types of driver 
gases, pure Helium (He) and a mixture of 70% Helium and 30% 
Nitrogen (N,), driven into Nitrogen as a test qas. An 
exploratory shock-tunnel nozzle experiment was conducted, 
and the capability of recording short-duration flow 
phenomena using available visualization techniques was 
demonstrated. 

It was concluded that the development of an operational 
hypersonic shock tunnel was feasible and would provide the 
means for follow-on efforts in the area of CFD code 
validation. At the same time, such a facility would provide 
a useful illustration of hypersonic flow to aeronautical and 
space applications engineering students. 

The following material is divided into three sections; 
1) a review of shock tube and shock tunnel ideas, 2) a 
discussion of the experimental investigation including a 
description of the facilities, experimental procedure and 
results; and (3) an evaluation of nozzle throat and dump 
chamber sizing requirements for a hypersonic tunnel. 


Instrumentation modifcations will be required in the 


development of a hypersonic capability. Appendix A reviews 
the existing facility and provides recommendations for 
upgrading the current hardware and instrumentation. 
Simplified flow models are developed in Appendices B апа с 
to enable the selection of nozzle throat diameter and dump 
chamber volume respectively. Experimental pressure and thin 
film traces obtained in the experimental program are 
contained in Appendix D together with predicted results 
obtained using a computer performance code, which is listed 


in Appendix E. 


паа ШФНЕОКЕТГСАР CONSIDERATIONS--A REVIEW 


A. BASIC SHOCK TUBE 

A shock tube is an economical short duration test 
facility used in laboratory experiments in gas dynamics, 
aerodynamics and physics. The principles of shock tube 
operation , which are well documented in numerous texts, for 
example, Gaydon and Hurle [Ref. 2] and Owczarek [Ref. 3], 
are illustrated in Figures 2 and 3 (Ref. 4]. As shown in 
Figure 2, the device consists of a long tube divided into a 
high pressure and a low pressure section by a diaphragm. 
When the diaphragm separating the high pressure driver 
section, region 4, and the low pressure section, region 1, 
is ruptured, a series of compression waves coalesce into a 
single shock front (S) which compresses and heats the region 
1 test gas to region 2 conditions. Simultaneously, a series 
of rarefaction (expansion) waves (R) travel into the high 
pressure section, region 4, cooling and expanding the driver 
gas to region 3 conditions. 

Region 3 (the expanded driver gas) and region 2 (the 
compressed test gas) are separated by a contact surface (C). 
The gases on either side of the contact surface travel with 
the same velocity, and are at the same pressure, but are 
different in temperature. A time history of the process 


(time vs. distance) is shown in Figure 3. 
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What is shown in Figure 3 can be determined accurately 
using the method of characteristics and shock theory (see 
[Ref. 3]). As the shock front reflects off the driven tube 
endwall, the region 2 gas is effectively shocked a second 
time, creating in region 5 stagnant gas at high pressure and 
temperature. This region may be used for high temperature 
gas dynamic studies or, when coupled to a nozzle, made to 
Serve as a reservoir of compressed and heated gas for a 


hypersonic wind tunnel. 


B. PERFORMANCE 
The strength of the generated shock (M) is related to 
the pressure ratio (P,/P,) across the diaphragm and the speed 
of sound in the undisturbed gases through the equation given 

by Gaydon [Ref. 2]. 
2 


2> 


-(- — 
Ave v Pon — X ===> “/ = a Y -] 
Р. (1+1) (7+1) ад ШЕ 


Flgures 4 [Ref. 5] and 5 [Ref. 6], which illustrate Equation 
(1), show that the strongest shocks are generated by using 
high pressure ratios and by using driver gases with high 
speeds of sound and low ratios of specific heats. 

As the desired shock strength increases, the required 
P,/P, increases significantly. To achieve stronger shocks 


several approaches are possible. The most common include 


Р /Р 
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Figure 4. Pressure Ratio (P,/P,) Required to 


Generate a Given Shock Strength Using 
Various Gas Mixtures [Ref. 5] 
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Figure 5. Variation of Shock Strength as a Function 
of Speed of Sound 
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(b) Driver Gas Specific Heat of 1.67 and a Driven 
Gas Specific Heat of 1.67 


Figure 5 (CONTINUED) 


TI 


the use of an area reduction between the driver and driven 
sections [Ref. 5], the use of multiple diaphragms with a 
buffer gas in the driver [Ref. 5], and heating the driver 
gas through combustion or electrical discharge techniques 
[Ref. 6]. Figures 6 and 7 [Ref. 7] show the relationship 
between shock strength (M.) and reflected region 5 


conditions for equilibrium air. 


Gi SHOCK TUNNELS 

The reflected region 5 pressure and temperature may be 
used to supply a supersonic or hypersonic nozzle if the 
nozzle throat is made sufficiently small to ensure full 
reflection of the incident shock. The evolution of this 
arrangement is shown schematically in Figure 8. To maximize 
test time and minimize the nozzle starting time, often a 
diaphragm separates the nozzle and dump chamber from the 
driven end of the shock tube. The dump chamber is sized and 
evacuated to a level which ensures that, when the tube is 
fired, the pressure in the chamber remains below the static 
pressure at the nozzle exit for the duration of steady flow. 
[Ref 8:pp. 624-640] 

For hypersonic applications a conical nozzle is often 
used. Table 1 provides calculated test conditions for 
various flow expansions using ideal gas theory. The 
computer program used for the calculations is described in 


Appendix E. Since, in an isentropic expansion, static 
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Reflected Pressure as a Function of Shock 


strength for Equilibrium Air [Ref. 
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Figure 7. Reflected Temperature as a Function of Shock 
Strength for Equilibrium Air [Ref. 7] 
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ТАВІЕ 1 


THEORETICAL PREDICTION OF TEST CONDITIONS FOR A HYPERSONIC 
NOZZLE (M, = 10) DRIVEN FROM THE REFLECTED REGION 
OF A SHOCK TUBE 


Driver Р. Т. Р, T: рашае 
Gas M. (Atm) (°K (mm Hg) (°K) (msec) 
2.4 58 ove 1.02 50 . 041 
37 32 2000 700 94 029 
Helium 4.2 23 2460 .41 116 42026 
4.5 19 2700 а IST .024 
26 23 632 . 4 30 JOST 
2.4 oM 91 . 169 50 . 041 
Nitrogen 
26 o 1110 - 116 52 2099 
265 5 3 1200 209 56 Оби 


temperature decreases as the Mach number (of the test 
sectlon) increases, liquefaction of the test gas becomes an 
issue. Since nitrogen closely matches the aerodynamic 
properties of air, but liquefies at a lower temperature, it 
is preferred for hypersonic facilities in which simulated 
real gas effects are not required. This is illustrated in 
Figure 9 which also shows the interesting results of wind 
tunnel experience [Ref. 9]. In practice, condensation is 
observed to occur at temperatures which are lower than the 


saturation temperature, thereby reducing somewhat the 
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stagnation temperature requirements to reach a given Mach 
number. 

Given the extremely short test time available in a shock 
tube, any coupling to a hypersonic nozzle must be done in 
such a way so as to minimize nozzle starting transients and 
provide the maximum available steady flow in the test sec- 
tion. The starting time, defined as the time to approach 
steady flow through a hypersonic conical nozzle (7, = 7/5), 


1s given by Smith [Ref. 8] as 


768 
SS. 0 2E E a 
torpe کے‎ ся м“ ІҢ «ам Ў) БМ.” ^-1)] (2) 


A reduction in starting time may be achieved by using a 
large nozzle angle (ñ = nozzle half angle), a small throat 
radius (Е.) and a high stagnation temperature in the 
reflected region.  [Ref. 8:pp. 625-640] The variation in 
starting time with increasing incident shock strength for a 
conical nozzle with ñ = 30 degrees is shown in Table 1. 
After the unsteady starting process, the useful test 
flow will continue as long as a sufficiently high pressure 
ratio is maintained across the nozzle. The steady flow 
process will be a function of only the upstream reflected 
stagnation conditions as long as the dump chamber has 
sufficient volume. Ideally, uniform reservoir stagnation 
conditions will end when the reflected region conditions are 


changed by a wave emanating from the interaction of the 


15 


contact surface and the head or tail ої the reflected expan- 
Sion wave (Figure 2). The useful test time in the nozzle 
ends when the reflected pressure drops below the level 
required to maintain the static pressure at the nozzle exit 
above the back pressure in the dump chamber. Thus the 
useful test time will depend on the nozzle throat size, 
which governs the flow rate out of the tube, and the size of 
the dump chamber, which determines the rate at which the 


back pressure builds up. 


D. TAILORED OPERATION 
D Tailoring 

The fundamental limitation of a shock tube is the 
avallable test time. Considerable effort has gone into 
trying to extend shock tube test times [Ref. 10]. Test 
times can be lengthened by as much as a factor of ten by 
"tailoring," which means adjusting conditions so that the 
ratio of the speed of sound and specific heat behind the 
reflected shock matches those in the region behind the 
contact surface. When this tailored condition is present, 
the reflected shock passes through the contact surface 
bringing it to rest. A Mach wave is reflected forwards off 
the contact surface and generates no additional disturbances 
as it passes through region 5. In the untailored condition, 
the wave reflected forward from the contact surface is not a 


Mach wave but a rarefaction wave when a, < a}, and a shock 


1. 


wave when a, > ay. The reflected waves cause conditions in 
region 5 to change. Experimental evidence of the effect of 
the wave interactions obtained by Pennelegion [Ref. 10] are 
shown in Figure 10. Pennelegion's research fully explored 
tailoring analytically and provided supporting experimental 
evidence. His findings and experimental techniques were 
followed in the experiments undertaken in the present study. 
The tailored interface technique requires that the following 


equation be satisfied [Ref.10:p. 4]; 


y, P. 
“2 _ 1. атаны (3) 
а, ү үші P 
ЕЕ 
"ТҮ? 


A unique tallored condition (shock Mach number), exists for 
every combination of driver and driven gas. The duration of 
constant pressure in region 5 lasts until the rarefaction 
wave reflects from the driver endplate and returns to meet 
the contact surface. For low incident shock strengths it is 
the head of the rarefaction wave, and for high incident 
shock strengths it is the tail of the rarefaction wave which 
terminates the constant pressure. 

The tailored interface technique for producing 
extended periods of constant reflected pressure and 
temperature was demonstrated experimentally for Mach numbers 


below 3.4. When the technique was applied to shocks of Mach 
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Figure 10. Experimental Observation of the Three Types 
of Wave Interactions in the Reflected Region 
of a Shock tube [Ref. 10] 
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number greater than 3.8, there was a significant drop in 
test time due to boundary layer interaction. In the under- 
tailored condition the reflected conditions showed an 
initial fall in temperature and pressure, followed by a slow 
rise in pressure and continued drop in temperature. The 
overtailored condition revealed a pressure rise from the 
reflection from the contact surface followed by a rapid drop 
in temperature. In all cases with strong incident shock 
waves, the real gas effects caused significant variation in 
the reflected region conditions.  [Ref. 10] 

2. Mixing of Driver Gases (He and N,) 

The application of tailoring to increase test time 
in a hypersonic shock tunnel allows more accurate 
measurements of force and pressure to be made on test 
models. As desired test section Mach number increases there 
1s a significant decrease in working section static 
pressure. It is therefore desirable to have reflected 
region pressures as steady and as high as possible. Pressure 
recovery, defined as the ratio of the initial step in the 
reflected pressure (Р.) to driven pressure (P,), may be 
increased by increasing the speed of sound of the driver if 
some test time is sacrificed. The combination of a driver 
of Helium (He) and Nitrogen (N;) can provide an improvement 
in test time and pressure recovery over a driver of К, 
alone. Experimental evidence of this is shown in Figures 11 


and 12 from Pennelegion [Ref. 10]. From the same study, 
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Figure 11. Available Steady Pressure Time [Ref. 10] 
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Figure 12. Experimental Pressure Recovery for Different 
Driver Gases [Ref. 10] 
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Figure 13 and Table 2 show tailoring using various mixtures 
of He and N, to generate shocks in the range 1 < M, < 3.4. 
The tailoring shock Mach number may be selected to provide 
the required stagnation temperature to prevent liquefaction 
and maximize test time for hypersonic nozzles with various 
Mach numbers.  Pennelgion (Ref. 10], provided the method for 


calculating tailored conditions for He, N, mixtures. 


4 я 5 Ф 
Mach number for tailoring 





Helium traction 


Figure 13. Variation of Tailored Mach Number with 
Helium Fraction (Ref. 10] 
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TABLE 2 


EXPERIMENTALLY TAILORED MACH NUMBERS FOR DRIVER 
GAS MIXTURES [REF. 10] 


He Time Steady T. 
Fraction Mos P,/P, Pressure (msec) (°K) 
ш. 7 290 24 635 714 
0.8 27 60 6.0 E05 
0.9 2.8 70 4.7 1159 
1.0 3.4 124 4.3 03 
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Im this section the existing shock tube and instrumenta- 
tion, and the modifications made during the course of the 


present experiments, are described. 


A. STRUCTURE 

The three inch diameter shock tube at the Naval 
Postgraduate School Aeropropulsion Laboratories is shown in 
Figure 14. A detailed description is given in Penaranda 
(Ref. 4]. The tube, constructed of highly polished 303 
stainless steel consists of a driver section 7.75 feet long, 
shown in detail in Figure 15, and a driven section 18.9 feet 


Tong. 


B DIAPHRAGM CHAMBER 

The tube incorporates a double diaphragm arrangement to 
maintain precise and repeatable diaphragm bursting pressure 
ratios (P,/P,). The tube currently is fitted with diaphragms 
designed to burst using a constant driver pressure of 600 
psi. Scored brass diaphragms designed to rupture at 400 psi 
are inserted. Variation in shock speed is obtained through 
adjustment of the downstream pressure (P,). To obtain the 
desired pressure ratio, the 2.5 inch long diaphragm chamber 
(Figures 15 and 16) and the driver chamber are simul- 


taneously filled to one half of the desired P4 pressure (300 
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Figure 15. Shock Tube Driver Section (Ref. 4] 
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Figure 16. Dual Diaphragm Chamber (Unbolted) 


psi). A diverting valve is closed to isolate the diaphragm 
chamber at 300 psi. The driver section pressure is then 
raised to the desired P, pressure (600 psi). To fire the 
tube, a diverting valve is opened and some of the driver gas 
at 600 psi is dumped into the diaphragm chamber causing the 
first diaphragm to break. The second diaphragm breaks im- 
mediately, as the pressure on the driven face of the 


diaphragm rapidly falls. 
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The dual diaphragm arrangement minimizes bursting 
pressure fluctuations due to inconsistencies in the material 
of the diaphragms. By applying a pressure of 200 psi above 
the design bursting pressure of the first diaphragm, rupture 
of the two diaphragms is ensured with a precise value of P,. 
In shock tubes using a single diaphragm arrangement, precise 
bursting pressures cannot be repeated unless the complexity 
of a diaphragm puncturing device is used. Without a 
puncture device, the bursting pressure fluctuations remain a 


function of the uniformity of the diaphragm material. 


C. NOZZLE ASSEMBLY 

An existing two dimensional stainless steel nozzle 
designed for an exit Mach number of 5.72 (A/A* = 43.7) anda 
rectangular dump chamber to which it bolted were available 
for the present study. The assembly is shown in Figure 17. 
For the present experiments two 1" diameter window ports 
were installed to allow optical measurement techniques to be 
exercised, and a model support was constructed to position a 
3/8" axisymmetric model on the nozzle centerline in line 
with the window, where the Mach number was about 4.3. The 
final assembly is shown in Figure 18. The dimensions of the 
nozzle are given in Figure 19 and the design of the window 
installation is shown in Figure 20. Views of the nozzle and 
model support are shown in Figure 21. Since the nozzle and 


dump tank existed (from a previous gas dynamic laser 
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Figure 18. Assembled Dump Chamber, Model Support 
and Nozzle 
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Figure 19. Dimensional Drawing of the Nozzle 
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(а) Model, Model Support and Nozzle 
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Figure 21. Views of the Nozzle and Model Support 
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investigation), location and sizing of the window as well as 
the maximum flow duration were limited by the existing 


hardware. 


D. INSTRUMENTATION 
1. Shock брееа 

Shock speed was determined by measuring the time for 
the shock to travel between thin film gauges located at 
positions A and B and from B to the pressure transducer at 
C (Figure 14). The locally manufactured thin fllm gauges 
shown in Flgure 22, consist of a 0.25" diameter pyrex год 
coated with a 1/32" wide thin film of platinum giving a cold 
resistance of 100 ohms. [Ref. 4:p. 23] As the region of 
hot gas immediately behind the shock front passes over the 
platinum, its temperature and resistance change and a 
voltage pulse is generated. In the present work the signals 
from the thin films were amplified and recorded for 
analysis. The thin films at location A also served as a 
trigger for the digital counter and shadowgraph spark 
source. 

Shock speed was determined in different ways. The 
pre-existing techniques involved using oscilloscope records 
of thin film and pressure transducer outputs, and also the 
output of a digital counter triggered on and off by the thin 
films at A to B respectively. In the present work, in runs 


1 through 49 shock speeds were determined from a Polaroid 
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Figure 22. Thin Film Instrumentation Assembly 
[Ref. 4] 
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photograph of the oscilloscope, such as is shown in Figure 
23, and from the digital counter. Both techniques were 
identical to those used in the late 1960's. 

In runs 50-85, the outputs of the three traces were 
recorded digitally on three of the available four channels 
of a Scientific Alanta SD380 Spectrum Analyzer. Data were 
recorded at .019 and .039 msec intervals depending upon the 
time scale selected. An example is shown in Figure 24. The 
shock speeds were determined during postprocessing of the 
data. 

2. Driver Tube Pressure 


A dial gauge with 10 psi pressure increments was 


used to determine when the constant 600 psi driver pressure 


was reached. A dead weight tester was used to calibrate the 
gauge. 
3. Driven Tube Pressure 


Two calibrated Wallace and Tiernan Bourdon tube 
gauges were used to determine the driven tube pressures. In 
the range from 0 to 50mm Hg, a gauge with increments of .2 
mm Hg was used. For higher driven pressures, a second gauge 
with a range from O to 35 psi in increments of .05 psi was 
used. 

4. Driven Tube Temperature 
A Type K (chromel-alumel) thermocouple, connected to 


a Daigger Scientific Digital Panel Meter using matched 


91 





Figure 23. 


Typical Polaroid Photograph of the 
Thin Film Response from Location B 
and the Output of the Original 
Pressure Transducer 
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extension wires, was mounted in the driven end of the tube. 
The digital readout was recorded prior to every run. 
5. Reflected Pressure and Test Duration 

A piezoelectric pressure transducer, Endevco model 
2501, a charge amplifier, Endevco model 2700, and a 
Tektronics 551 oscilloscope were initially used to record 
the transient region 5 pressure. The lack of resolution and 
strange qualitative behavior observed in the results (Figure 
23), led to the mounting of an improved transducer. 

A PCB Piezontonics, Inc. Model 124A pressure trans- 
ducer connected to a model 504A dual gain charge amplifier 
were used in runs 50 through 75. The system was calibrated 
using a standard static calibration feature described in the 
manufacturer's manuals (Refs. 11,12]. The transducer 
provided a response of 240 psi/volts. Experimental traces 
obtained are given in Appendix D. 

6. Oscilloscopes 

Two Tektronics oscilloscopes, models 549 and 551, 
fitted with Polaroid cameras were used for runs 1 through 
49. The capability of the scopes are fully described in the 
associated manuals (Refs. 13,14]. Time measurements to 
determine shock speed, and voltage levels to determine 
pressure, were determined by overlaying a scale on the 


Polaroid photographs. 
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7. Signal Analyzer 


A Scientific Alanta SD380 four channel signal 
analyzer (Figure 25) equipped with a hard disk and plotter 
was used to record and process data from runs 50 through 85. 
The independent channel trigger delays (+), adjustable 
scales and post-processing capabilities of the digitized 
data (in contrast to photographing oscilloscopes), is 
evidence of the advances in measurement capabilities over 
the past 20 years. A full description of capabilities of 


the analyzer is given in the instrument manual. [Ref. 15] 


E. OPTICAL MEASUREMENT SYSTEMS 

A simple shadowgraph arrangement was generated from the 
components of a 6" schlieren system. The optical layout is 
shown in Figure 26 [Ref. 16]. The spark source was used 
with Polaroid film in a camera back attached to the nozzle 
window. The continuous source was used with a Hycam movie 
camera. 

1. Light Sources 

a. High Pressure Mercury Arc Lamp 
A high pressure mercury arc lamp (Per Type 112- 

5145) constant light source was used with a six inch 
achromatic convex lens with a focal length of 60 inches. The 
arrangement, mounted on an adjustable bench (Figure 27), 


provided a parallel light beam through the nozzle windows. 


41 


27 


Ж 
“Z 
А 


“%% 


[| 


2 
т; 
3 


Ж 


% 


сМ 


Figure 25. 


О ЕР 


оо УФ ИЯ А 


У АРА 


њи 


2 


E 
cate,’ che kd! 


г 2 


% 
ж E РРР, ж 


% 2 
UE 2 


7 و‎ 
Lm и 


4 £ 7 
ПОД СОЛО си 


2: 


L 4%; Зан a 


Н 
4 
> 
(4 


2% 


ПО НО СРС С 





E. 3 ر‎ 
et UU 
2 E 


Scientific Alanta SD380 Four 
Channel Signal Analyzer with 


Hard Disk Installed 


42 


! "| 0 


€— DRIVER SECTION 


&—- DIAPHRAGM CHAMBER 


€— DRIVEN SECTION 


6" CONDENSING LENS SHADOWGRAPH FOCAL PLANE 


- POLAROID FILM PLATE 
ре з LOCATION 


# ж=—===-—-| = с | ЕЕ НУСАН MOVIE 
| САМЕКА 
50" “яш 


€— DUMP CHAMBER 


- SPARK LIGHT SOURCE 
- MERCURY ARC LAMP HN 


Figure 26. Optical Layout Sketch 
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Figure 27. Constant Light Source and Lens Mounted 
on an Adjustable Optical Bench 
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Figure 28. Spark Source with High Voltage 
Power Supply 
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The beam was adjusted by the procedures described in the 
assoclated manual [Ref. 16]. 
b. NASA 229-D Spark Light Source 
A NASA 229-D spark light source (Figure 28) 
provided a nanosecond spark for shadowgraph pictures taken 
with Polaroid film. The spark was triggered from a thin 
film connected through the delay and voltage circuits shown 
in Figure 29. The timing for the spark source was 
determined experimentally. 
2. Camera Systems 
a. Hycam Movie Camera 
A Hycam movie camera [Ref. 17] with a framing 
rate on the order of 6000 pictures per second at runout of a 
400 ft roll of film, was used. The camera was equipped with 
a Nikon 55mm lens and ASA 250 black and white film. A time 
reference наз provided on the edge of the film by a strobe 
light pulsing at 1000 pulses per second. The camera was 
triggered manually. The camera was mounted on a second 
adjustable bench as shown in Figure 30. 
b. Spark Shadowgraphs 
Spark shadowgraphs were taken using a Polaroid 
Land Film Holder with Type 52 ASA 400 Black and White filn, 
held directly behind the optical windows. To eliminate 
stray light, the mounting arrangement was contained ina 


cardboard box, and tests were conducted after dark. 
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Figure 29. Spark Timing Delay Circuit and High 
Voltage Pulse Generator 





Figure 30. Hycam Movie Camera with Separate 
Pulse Generator 
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ТУ. EXPERIMENTAL PROCEDURE 


A. PREPARATION 

Prior to testing, with the driven endplate removed and 
the diaphragm chamber opened, a cloth lightly dampened with 
Acetone was pulled through the driver and driven sections. 
O-ring grooves of the opened sections where thoroughly 
cleaned and a thin coating of vacuum grease was applied to 
the rubber O-rings. Diaphragms were wiped clean of visible 
dirt and the system reassembled. The cleaning procedure was 
repeated after approximately 15 runs. Using a roughing 
pump, a clean-tube vacuum pressure of 3 mm Hg was routinely 


achieved. 


B. LOADING THE DRIVEN GAS 

Prior to loading to the desired test pressure (P,), the 
driven section was purged. To purge, the tube was evacuated 
to 3 mm Hg and then isolated from the roughing pump. 
Nitrogen was introduced to a pressure 1 to 4 psi. The tube 
was then evacuated a second time and repressurized to the 


desired test condition. 


C. LOADING THE DRIVER GAS AND FIRING 
The driver assembly is designed to operate at 600 psi 
with single or dual gas combinations. To purge the driver 


section and diaphragm chamber, the system was pressurized to 
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approximately 30 psi and then vented to atmospheric pres- 
sure. The system was then pressurized to 600 psi with the 
diverter valve to the diaphragm chamber being closed at 300 
psi. When the 70% He and 30% N, driver mixture was used, 
the system was pressurized to 420 psi with He and then to 
600 psi with N;. 

The tube was fired by opening the diverter valve to the 
diaphragm chamber. Some seepage from the driver section was 
evident. To minimize this effect the tube was fired 
immediately after reaching the 600 psi level. Replacement 
of O-rings would minimize this problem.  Seepage into the 


driven section was not evident. 
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У. EXPERIMENTAL RESULTS AND DISCUSSION 


A series of experiments was conducted to gain experience 
with the shock tube. The study consisted of approximately 
85 shots involving the use of He and a 70% He/30% N, mixture 
as driver gases, driven into a N, test gas. Studies with 
the shock tube assembled with the pressure transducer 
mounted in the end of the driven tube, were concentrated on 
investigating the concept of tailoring to maximize the 
magnitude and duration of steady pressure in the reflected 
region. While more than 75 shots were made in this 
configuration, the replacement of the pressure transducer 
and postprocessing of data taken with the SD380 greatly 
enhanced the results. For this reason, only data taken after 
run 50 using the SD380 are presented. Studies with the 
nozzle attached at the end of the driven tube were 
exploratory attempts to obtain flow visualization in a shock 
tunnel configuration. The results of the two series of 


experiments are discussed in the following sections. 


A. TAILORING 

Using the nomenclature of Glass [Ref. 18], a tailored 
condition is defined as the condition when a Mach wave is 
reflected off the contact surface and causes no change in 


the reflected pressure (P.,). This can be identified 
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experimentally (Figure 10) as the condition when the ratio 
of the initial reflected pressure (P,) and the pressure 
after the wave interaction at the contact surface reflects 
off the endplate (P,) being equal to one. Figure 31 shows 
examples of results in the present study for the end wall 
pressure untailored and tailored conditions using He with an 
incident shock strength, M, = 3.4. Similar traces for near 
to a tailored mixture (.7 He and .3 N;) with M, = 2.0 are 
shown in Figure 32. These results are in complete agreement 


with Pennelegion [Ref. 10]. 


B. TIME OF STEADY PRESSURE 

It was experimentally determined in the present study 
that a tailored mixture as a driver gas produced a longer 
duration of steady pressure than that of a pure He driver. 
The time of steady pressure for the existing NPS shock tube 
using He as a driver was approximately 3.6 milliseconds 
while the mixture extended the test time duration to 5 
milliseconds. А comparisons of the test times obtained in 
the present study with those of Pennelegion (Ref. 10] is 
shown in Figure 33. The trends are seen to be similar. 
Since the work of Pennelegion (Ref. 10] was performed at a 
constant driver pressure of 1000 psi in a N.P.L. 2" shock 
tube with a 10 ft driver and a 12 foot driven section some 
differences in the magnitudes of the test times are to be 


expected. 
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Figure 31. Experimentally Determined Endwall Pressure 
Traces Using Helium as a Driver Gas 
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Comparison of Steady Pressure Time Achieved 
in the NPS 3 Inch Shock Tube with Experi- 
mental Results of Pennelegion (Ref. 10] 
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For low strength shocks the head of the expansion wave 
causes an early degradation of the reflected pressure in the 
NPS tube. The longer drlver section in the N.P.L. tube is 


responsible for the difference. 


C. MPRESSURE ВЕЕТ VERY 

Figure 34 shows the fraction of initial driver pressure 
available in the reflected region as a function of shock 
Mach number. The data show that as shock strength increases 
pressure recovery decreases for both He and mixed driver 
gases. The data are in good agreement with the results of 
Pennelegion (Ref. 10]. The measured values of the reflected 
shock pressure, which are the pressures available to drive a 


hypersonic tunnel, are given in Tables 3 and 4. 


D. “REFLECTED TEMPERATURE 

The reflected temperature was not determined 
experimentally. Instead, shown in Table 2 are values of 
stagnation temperatures calculated for shock strengths in 
the range from 2.0 to 3.4. These values were taken from the 
tabulated solutions (based on equilibrium chemistry) of 


Bernstein [Ref. 19:рр. 32-43]. 


E. SHOCK ATTENUATION 
A representative illustration of measured shock 
attenuation for both a pure gas driver and the driver gas 


mixture is shown in Figure 35. Velocity can be inferred 
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Figure 34. Comparison of Reflected Pressure Measurements 
іп the NPS 3 Inch Shock Tube with Experimen- 
tal Results of Pennelegion (Ref. 10] 
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TABLE 3 


EXPERIMENTALLY DETERMINED REFLECTED PRESSURE MAGNITUDE 
AND DURATION AND CALCULATED TEMPERATURE (HE DRIVER/N, DRIVEN) 


Run Р. Time of Steady 
# М, BAP; (psi) Pg/ P5 P, (msec) 
50 3.08 994 150 go 123 
Ба 2592 353 235 1.25 1.2 
52 4.39 414 207 2-39 19 
53 3.32 93 ЧЕ 200۳87 125 
Dd 4.52 500 183 125 0 
58 защо 207 304 ш. 23 О 
59 3.66 148 386 1.13 1.6 
60 3:36 169 264 1213 З то 
61 972 167 317 ЛЗ 0:25 
64 3.42 164 == 
65 35.65 155 == 
66 184 
67 dE ا‎ 
68 291 400 
74 4.50 667 
75 4.10% 1240 


Estimated Value‏ ٭ 


TABLE 4 


EXPERIMENTALLY DETERMINED PRESSURE MAGNITUDE AND DURATION 
AND CALCULATED TEMPERATURE (70% HE, 30% N, DRIVER/N, DRIVEN) 


Run Р. Time of Steady 
i M, PUE (psi) Pg/ P5 P, (msec) 
62 2.40 50 3:17 1.10 2.8 
63 2.66 51 310 1523 l9 
69 2.47 50 310 1.09 2:85 
70 2525 50 305 quon 3.3 
71 29803 24 427 252 4.9 
72 2.03 22 422 1.08 5 
73 1.99 21 410 ies 2 
77 
78 155 
79 3.0% 171 
80 3.4% 638 


Estimated Value‏ ٭ 
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from the slope of the lines between the thin films А апа В 
and between B and the pressure transducer at C. An 
increasing slope indicates a slower speed. The general 
trend is that for a given gas mixture, the stronger the 
shock strength the smaller the attenuation between B and C. 
Reflected velocity is of the order of one half the incident 
velocity. Tabulated results for all runs are found in 


Appendix D. 


F. OPTICAL MEASUREMENTS 

Before attempting to obtain photographs during the 
firing of the shock tube, a small freejet from a 1/4" 
flattened copper tube, powered by shop air, was installed in 
the model location, the spark light was adjusted using a 
neutral filter and a successful shadowgraph was taken. 

1. Shadowgraph Using Hycam Camera 

The Hycam Movie camera was triggered manually either 

one or two seconds before the diaphragm was ruptured. The 
400 ft rolls of ASA 250 speed film provided approximately 
one second for the camera to accelerate, and a three second 
window to capture the event. It was found that all the 
available light was needed to obtain a correct exposure and 
the camera was focused on the outside surface of the window. 
The shadowgraph was made visible by the faint film which by 
then had deposited on the inside of the window. The result 


was successful. 
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Wlth the blunt body shown in Figure 21(a) installed, 
the nozzle starting Mach wave (Figure 36(a) and (b)), and 
the establishment of steady flow with an oblique shock 
visible (Figure 36(c)), were successfully recorded. Steady 
flow (started) conditions lasted a brief time until dump 
chamber back pressure unstarted the nozzle. Incident shock 
strength was estimated to be M, = 3.8 providing a peak 
pressure in the reflected region of approximately 460 psi. 
The actual shock speed was not measured during the test 
Since the starting of the camera caused premature triggering 
of the spectrum analyzer. However the initial conditions, 
and therefore the reflected pressure, were the same as in 
run 58. A complete record of this run given in Appendix D. 

2. Spark Shadowgraph 

Before firing the shock tube, the correct exposure 
for the expected shadowgraph was verified by opening the 
film and firing the spark source. The result is shown in 
Figure 37. It was subsequently found that when the shock 
tube was fired, the thin film signal failed to trigger the 
spark, and thus any exposure on the film was the result of 
light produced from within the tube. The results of 
exposing the film without discharging the spark as the shock 
Mach number was progressively increased, are shown in 
sections (a) to (e) of Figure 38. A faint glow appeared 
first аға М, - 3.6. As M, was increased for both a pure He 


driver and mixed driver gases, the luminosity increased. 
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Figure 36. Results Obtained Using Hycam Camera 
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Figure 37. Spark Shadowgraph with No Flow Present 
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The source of the light requires further investigation, but 
1t is clear that some contamination must be present since 
the driver and driven gases, if pure, would not be excited at 


the temperature produced in the firing process. 
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УТ. SIZING CONSIDERATIONS FOR A HYPERSONIC TUNNEL 


The experimental program showed that the basic shock 
tube gave reliable and repeatable results and, following the 
upgrading recommendations detailed in Appendix A, it could 
well serve as a driver for a hypersonic tunnel. Having 
experimentally defined the stagnation temperature and 
pressure available to drive such a facility, it only 
remained to arrive at a nozzle and a dump chamber 
configuration suitably matched to the tube. Simplified 
models were therefore developed to, first, calculate the 
effect of nozzle throat diameter on the available steady 
flow time and, second, to determine a suitable dump chamber 
Size. The two calculations are detailed in Appendix B and 


Appendix C respectively. 


A. NOZZLE THROAT SIZE 

The nozzle throat diameter must be sized such that the 
time for which conditions in the reflected region are steady 
must be greater than the time required for the compressed 
mass of gas to exhaust through the nozzle. The model 
described in Appendix B estimates the time in relation to 


the falling pressure in the reflected region to be given by 


тн zx (B.10) 


с "ATP: Y сз ee 


If it is assumed the supply pressure can drop to half its 
initial value, the time required can be calculated using 
Equation (B.10) for different driver gases and shock 
strengths. Through the ratio по/ то, the calculated time is 
implicitly a function of the throat diameter. Results of 
such calculations are shown in Figure 39(a) and (b). The 
four curves shown illustrate the dependence of useful time 
on the shock strength and throat diameter. Clearly, the 
shock strength is less important than the throat area. 
Based on these results, a throat diameter 6.5mm (.025 
inches) was selected. This choice made best use of the 
steady flow time which was experimentally determined to be 
from 3.5 to 5 milliseconds using tailored conditions. For 
the 6.5mm throat, Equation (B.10) was used to calculate the 
time for the reflected shock region pressure to decay to 
various fractions of its initial magnitude. The results are 
shown in Figure 40. Clearly, for even the strongest shock, 
the time for the pressure to decay to 80% of its initial 
value is longer than its experimentally determined tailored 


Fest time of 5 msecs. 


Eee DUMP CHAMBER SIZE 
The dump chamber requires a capacity such that the back 
pressure remains below the nozzle exit pressure for the 


duration of steady conditions in the driven tube. The 
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analysis in Appenix C gave the required ratio of dump 


chamber volume to reflected shock reservoir volume as 


ре ат. 
0 Банні 0 5 s 
P 





Eo D. 
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Un 


Examination reveals that the parameters which most influence 
the required dump chamber volume are the incident shock 
strength (which determines the reflected reservoir 
conditions), the desired test section Mach number, the 
driver gas selected, and the specified reflected region 
pressure (B,/B,j) at which the flow breaks down. Since the 
pressure recovery is a function of which driver gas is 
selected, calculations were made for both nitrogen and 
helium. The results are shown in Figure 41. A mixture of 
the driver gases would be expected to fall between these two 
curves. Consistent with the finding above, the criterion 
that the dump chamber pressure becomes equal to the nozzle 
exit pressure when 20% of the stagnation pressure was 
exhausted, was taken to determine the necessary dump chamber 
volume. Shown in the figure are results obtained assuming 
the unknown temperature in the dump chamber to have its 
extreme maximum and minimum values; namely, initial 
reflected shock temperature and initial dump chamber 


temperature. From an examination of the results, a 
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Figure 41. 


conservative selection would be 15 cubic meters. This is a 


vessel approximately 6 feet in diameter and 15 feet long. 


C. HIGH ENTHALPY MODIFICATION 

Higher enthalpy flows can be generated using the free- 
piston arrangement shown in Figure 1. Stalker [Refs. 20,21] 
has carried out extensive experimental work in this area. A 
preliminary analysis shows that the NPS shock tube has the 
same internal diameter as Stalker's lst generation free- 
piston tunnel, but has a thinner wall by 1/2". Structural 
limitations would prevent the full duplication of Stalker's 
tunnel but, with the published information on the design of 
free piston drivers, it should be possible to extend the 
present tube to an increased enthalpy level at which real 


gas effects of interest can be studied. 
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VII. CONCLUSIONS AND RECOMMENDATIONS 


The present study showed that: 


- The existing shock-tube is straightforward and 
relatively inexpensive to operate, and gives 
repeatable test conditions; 


- The existing shock tube can be used to drive a 
hypersonic shock tunnel to a Mach number of 12; 


- Tailored operation can provide uniform stagnation 
conditions for periods of approximately 3.4 msec 
using Helium, and approximately 5 msec using a .7 
He/.3N, mixture, as the driver gas; 

- It is feasible to obtain both high speed movies and 
short duration spark flow visualization records of 
the test section flow field using currently available 
equipment; 

- Contamination in the tube can lead to unwanted 
radiation (which must clearly be eliminated if real 
gas effects are the purpose of the experiment). 

It is recommended that: 


- The existing tube be serviced and its operating 
instrumentation be upgraded, as outlined in Appendix A; 


- A Mach 10 nozzle, 15 cubic meters dump chamber and 
model support be designed and procured for operation 
with the current tube; 


- Consideration be given to the addition of a free- 
piston tunnel to provide high enthalpy simulation. 


Operation of a shock tunnel with the existing tube would 
provide a convenient and inexpensive facility for use in 
laboratory instruction. It would also provide a tool for 


thesis studies aimed at hypersonic CFD code validation. 
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Extension to a free-piston driver would provide an advanced 


tool for hypersonic research. 
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APPENDIX А 


APPRAISAL OF THES PRESENT SHOCK LUBE 


The NPS shock tube has been in operation since the late 
1960's. This appendix provides a description of recommended 
improvements/enhancements of the existing system, which 
result from the experience gained in the reported program of 


tests. 


А. О=к бо 

The planned relocation of the facility to the new Gas 
Dynamics Laboratory requires complete disassembly of the 
shock tube. Many O-rings are old and worn and they should 
be replaced at that time. Table 5 details the O-ring 


requirements. 


B. THIN FILM GAUGES 

The unavailability of replacement thin film gauges used 
to measure shock speed is unacceptable. Only the two in- 
stalled gauges remain operational. The options are either 
to manufacture thin films in-house as detailed in Penaranda 
[Ref. 4] or to procure them commercially through the supply 


system. 
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TABLE 5 
O-RING REPLACEMENT GUIDE 


Number Width of Depth of Inside Outside 
Location Required Groove Groove Diameter Diameter 








(inches) (inches) (inches) (inches) 

















Driver Tube 1 ‚ 100 . 8207 3.502+.005 
. 106 „163 > 
брасе/ 2 · 100 157 4.000 
Diaphraqm 106 М 4.005 
2 . 100 2154 3.500 
106 „163 Ваау 
Valve Body 1 101 AUT 3. 2485 
105 163 3. 2495 
Tubing 3 101 ОҒЫ? 3.2485 
Sections 105 105 3. 2495 
Instrumentation 5 101 E 3. 2485 
Section 105 ES 3. 2495 
Driver End Cap 1 101 157 3.2470 
105 163 3.2460 
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С. PRESSURE TRANSDUCERS 
1. Piezoelectric Transducer 

The two Endevco Model 2501 transducers with Model 
2642 M35 charge amplifiers should be returned to the manu- 
facturer for calibration. Biased and unbiased responses 
should be checked. Unknown, but apparently unrealistic 
response characteristics made these transducers unusable for 
the reported experiments. 

2. Piezoresistive Pressure Transducer 

Procurement of a piezoresistive pressure measurement 
system, including a transducer and amplifier such as the 
Kistler Model 603B Series and Model 5004 Charge Amplifier, 
would provided continuous measurements of shock tube 
pressures. 

A second transducer, to be used as a pitot probe for 
flow survey work in shock tunnel applications, is also 
required. Selection would be based on the selected 
nozzle Mach number and probe geometry. Cornell Aeronautical 


Laboratory publications by MacArther should be consulted. 


D. THIN FILM SIGNAL AMPLIFIER 

The current thin film signal amplifier has no record of 
calibration for the past ten years. Voltage fluctuations 
are apparent and dial adjustment response is very poor. 
Prior to the next use the amplifier components should be 


tested and replaced as required. 
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Eee DATA ACQUISITION 

It is recommended that the Scientific Alanta signal 
analyzer or its equivalent be used in all future testing. 
Post processing capability, and the availability of four 
channels with independent trigger adjustment provide a 
valuable experimental capability not available with the 


Tektronics 549 and 551 oscilloscopes. 


E. DIAPHRAGMS 

The existing brass diaphragms allow driver operation at 
a constant 600 psi level. Redesign to a higher operational 
level of approximately 700 psi would increase the range of 
reflected pressures available to drive a shock tunnel, with 
an associated cost of higher gas consumption. While not 
essential, this modification would reduce problems which 
would be expected in an extension to a shock tunnel with 
respect to the measurement of small aerodynamic forces. 

пол (Ref. 2:p. 91} and Glass [Ref. 18: pp. 503-511] 
provide valuable design insight into successful diaphragm 


material selection and preparation. 


G. MISCELLANEOUS HARDWARE 

Valve 7 and valve 11 located in the driven end of the 
tube [Ref. 4:Figure 2] exhibited some signs of leakage 
during the loading process. Repair or replacement of the 


valve seats is recommended. 
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APPENDIX B 


ESTIMATION OF SHOCK TUNNEL NOZZLE THROAT AREA 


In order to size a suitable nozzle to be driven by the 
high pressure and temperature in the reflected region, a 
Simplified model was developed. The model was based on the 
following assumptions: 


- A constant volume reservoir corresponding to the in- 
itial reflected shock conditions was available to drive 
the nozzle flow. The reservoir volume was defined by 
the diameter of the tube and the axial distance from 
the endplate of the tube to the point where the contact 
surface and reflected shock met (Fig. 3); 


- The gas remaining in the reflected region expands 
isentropically to fill the constant volume; 


- The flow through the nozzle throat was isentropic and 
always choked. (This implies that the evacuated dump 
chamber had sufficient volume, so as not to influence 
the nozzle flow.) 


The total mass initially in the reflected reservoir is 








given by, 
Бе 
m. = М» (Ва у) 
0 в 0 
where 
1D S 
М = С 
0 4 


The mass at any later time is given by 
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m ЕНН 


ш = SV (B.2) 





As the mass flows through the choked nozzle the rate of 


change of mass in the volume is equal to that which flows 


out through the nozzle, thus 











me 
2 (571) 
A*P Y ۷ء۰٢‎ 
Ee - оү where Су = в 2. (Bes 3) 
VIN g 
L 
Integrating 
Е t АТР, 
0 0 ШЕ 


Since the expansion in the reservior is assumed to be isent- 
ropic, the stagnation pressure and temperature are related 


a cording to 








P P 
ПЕ ғаны (В.5) 
Ys Ys 
ШЕ С 


Taking the ratio of Equations B.1 and B.2 and substituting 


into Equation B.5, 
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т 555 Р 1/ У 5 
мо (4%) Щи ше (B.6) 
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Substitution into Equation B.3 glves 








Ygti 
P Y 
. dmn  , ,)ٰ٤ ٦٣ 5 (ВЕНА 
е упо Lo He 
5 
so that 
m М 
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On integrating, the variation of mass in the reflected 


reservior with time is given by 





үсті 
m eL m : 
0 5 2 0 
те 


and the variation of pressure in the reflected region with 
time is given by 
П 
жй, майра О 9 I tL A (B.10) 


А 5 
С.А Р. Үс 1 
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To estimate a suitable nozzle throat, Equation B.10 was 
programmed and a ratio of P,/P, equal to one half was 
selected as the criterion for calculating the useful run 
time. Using throat area as a variable all other parameters 
were a function only of the incident shock strength and of 
the physical size of the tube. Thus Equation B.10 enabled 
the available time for steady flow through a given throat 


diameter to be estimated. 
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APPENDIX C 


ESTIMATION OF THE REQUIRED DUMP CHAMBER VOLUME 


In a shock tunnel application, the dump chamber pressure 
is required to remain below the nozzle exit static pressure 
for the useful duration of the reflected-shock supply condi- 
tions in the shock-tube. То aid in sizing a dump chamber 
for a Hava1 Postgraduate School shock tunnel the following 
analysis, based on a simplified model, was developed. The 
model assumes the pressure and temperature are uniform in 
the dump chamber, and that the velocity there is effectively 
zero. It is assumed that the flow to the exit of the nozzle 
is isentropic, and that the volume of the nozzle is negligi- 
ble in comparison to that of the chamber. The volume of the 
dump chamber is given in terms of the initial mass contained 


and its properties, ру 


D (CS 





Dividing by the reflected reservoir volume, similarly writ- 
ten in terms of the initial mass contained and its proper- 
ties, the ratio of dump chamber to reflected shock volume is 


given by 


р. р з В (ç 





As the flow continues, the original reflected pressure (P) 
decreases to lower values (P,). Since the nozzle area ratio 
is fixed, the pressure at the exit of the nozzle (P,) also 
steadily decreases. Meanwhile, the pressure in the dump 
chamber (P) steadily increases until, at a point in time, 
P, = P., and the useful test time is over. The ratio of the 


initial reflected pressure to dump chamber pressure is then 





given by 
25 = E = G5 gt (en 
ED Fe D t e 
where 
c 
Р үсті ED 
Баш 5 2 СА 
P o^ E 5 Ме ) 
e 
A mass balance gives 
ЕЕ z (G2 Ss) 
по mo m. T Myo 
so that 
m I. m 
ED e pa З. DO (о 
m. mo т 


Using Equation B.6 from Appendix B 
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Combining Equations C.2 through C.9 


(C 10) 





Equation C.10 provides a relationship between Vyj/V, and the 
value of P,/P, which terminates the useful test time, in 
terms of given initial reflected reservoir and dump chamber 
conditions. Using an interval halving technique the dump 
chamber volume requirements were investigated using an 
assumption that the temperature in the dump chamber was 
equal to its initial value, and also with the worst-case 
assumption that it was equal to T,. (It is noted that the 
temperature change can be included explicitly through the 


conservation of energy. However the complexity of doing so 
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15 not warranted when an estimate of the needed volume is 


all that is requlred.) 
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APPENDIX D 


SUMMARY OF EXPERIMENTAL DATA 


This appendix contains data taken using thin film gauges 
and the pressure transducer. First, in Figure 42, the 
measurements obtained of shock Mach number using oscillo- 
scopes (preceeding Run 50) are shown. For all runs from Run 
50 on, data were taken using the Scientific Atlanta SD380 
signal analyzer. Tables 6 and 7 list data for runs using a 
mixture of helium and nitrogen and for runs using pure 
helium as a driver, respectively. For each run, trace A 
recorded the response of the thin film at location A, trace 
B recorded the response of the thin film response at station 
B, trace C recorded the response of the endplate pressure 
transducer. Delay times were programmed prior to each run. 
Traces B and C were delayed by the same time. Figure 43 
provides definitions of items which appear on the SD380 
screen and identifies the menu selections of interest. 

The last 35 pages of the appendix contain printouts of 
the SD380 screen for the runs listed in Tables 6 and 7. 

Each printout can be interpreted with reference to Figure 
43. Times were obtained for shock travel between trans- 
ducers on traces A, B and C, using the time scale given on 
the printout, and the delay times listed in Tables 6 and 7. 


These results are included in the two tables. 
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APPENDIX E 


COMPUTER PROGRAM FOR SHOCK TUNNEL EVALUATION 


The attached program was used for the evaluation of the 
NPS shock tube and shock tunnel. Sources of shock tube 
relationships are noted in the program were applicable. 
Output was tabulated by the program and graphical analysis 
was provided by the plotting routine 'EASYPLOT' installed on 


the NPS IBM computer. 
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EDLE: 


CC 
pr 
СС 


СОСО СУ О 


رہ رق ری یہ ری ری С СУСС СУ О‏ 


DIDID 


SHOCKZ FORTRAN Al 


THIS PROGRAM CALCULATES THE IDEAL PERFORMANCE OF A SHOCK 

TUBE USING THE ANALYSIS PRESENTED IN THE ‘FUNDAMENTALS 

85۷۰۰ ت‎ ۱۱٦۷۷۱۱۸٦٦٦۷۱5 ۲۹۰۱۰٠ СТОИ ФАКЕК РО 595-598 (FIG 8-22) 

& ‘THE SHOCK TUBE IN HIGH TEMPERATURE CHEMICAL PHYSICS' BY GRA 
GRAYDON HURLE 1963. 


IMPLICIT REAL xX 4 (A-H, M, P-Z) 

REAL FCH XI 20, ХТОГ,ҒЕТОІ, ЕСНПІ,ЕСИ2 FONS, FCNG 

INTEGER ISNIN; J,;I1,I1I2 

EXTERNAL FCH,FCH1, FCN2, FCN3,FCNS 

DIHHEHOSIOIIFDTC200), oO (200);5C200),PGP10(200);POP10200) 
ИЛЛЕНЬТФПРЕЭРЯСОРО),М5(2002,05(2007, Т222000, Р5Р2(200),Р5Р6(200) 
РЕНО С5 (200), С6(20002Рорз200), 502 Шб/ТСС(200),ХС(200) 
ВОВЕ ШМЫУСС2(0О,/РӘРТСООО>ТЕТТ СОФОЛУТЕСОВОО,ВЕГТІНС200) 
DIMENSION ВТ 5102002, ТСС, ^_ 11200), DEUDAC200),PDICS) 

ЈУ АНОМ oc OO ХЕРЕ хръггс гориш хеогасгод)», хкбек! (200) 
DIMENSION ХТ2Т1( 200), 0605 (200), Т5ЭТАЕТС( 200), УРУОС8), МОС8), ААС8) 
ШГПЕПБЬТОПИРЕС ООО), ТӘСОИУЛТИКСОзОО),ШГСОШФЛУИТС200), ВТЕ5Т( 200) 
ШИПЕПЫТОЦ ET 0(200),110C200);1(C200, 6), TiASSCo00), DIDTHR(5),P1(0200 
ПАНЕ МО ТОМИ о 100,РІРТОСІЮ) 

60۲1110٣٣۸67 ۲٦ 018ء68 ۱۱۱۹۶ء‎ 

СОНМОН /А/ СОНУТ, САММА 1, САММА4 ,АЧА154, СЪ, С159, СС5, СС4, С5, АЧА1,С9 
60٦1۲10۱1٣ 67۹/۸۹۱۸5۱۴۰ ۰17 

COMMON / D7 ۳5۰ 

COMMON 7F7 PTPIOC,FT 

DATO Ol, FTOL, TF HEIN7Z0.1,6.0,0,0001,0.00001,0,507 


) 


— — «кн «ы, «шт «а: «а: «аз == == <= <= == <= = == <= = = = = = = = = <= = = = <= = = = = = = = = = = = = = = = = = = == = = = = = = = = = = 
- -. ee 


SELECT GAS CONSTANTS 


GAS MH GAMMA 5۲ OF SOUND 
AIR 28897 1.4 345 M/5 

ARGON 39794 167 1 

НЕ 4 1.67 1062 

12 28.02 1.4 250 

coz 44.01 1.229 269 


DRIVEN GAS CONSTANTS 
К1:8516/28.02 

11=288 

САММА1=1. 9 

۲12-69 


DRIVER GAS COHSTAHIS 
Ка-8514/28.02 

14-288 

САММА4=1. 90 
PG=600%6.89498E03 

A IS INCREMENT ON P4P1 LOOP 


А-100 

MM 9۹۲) OF LOOFS 
ММ=6 

РаР1(1)=20 


= = = = = = = = = = == = <= == = = <= = = = = = = = дю» = = = ле о жа “ша “ша ош» ee eee <= = = = = = = = = = = = = = = = = = = == = 
““ -н к» -- «ше «е «т «т «ші «ша «ж «ае «ша «ае «аз «ж» eee See ee ee 


INPUT TUNNEL PARAMETERSCDRIVER LENGTH L1 & L4) 
(1:6.455 
Eg 22505 
DT1=3%.0254 
TO EVALUATE TEST CONDITIONS SET FL2=1 
THR1=.0065 
FL1=0 
5ھ‎ CALCULATES TEST SECTION MACH # 


1 
й. IS FOR HOZZLE SIZIHG CALCULATIOHS 
Lo 
CALCULATES TEST SECTION CONDITIONS 
FLG=1 
1728 


FILE: 


O OO 


OOOO 


ہے" ہہ ہے та‏ - 


SHOCK2 FORTRAN Al 


FL6 EXPERIMENT OUTPUT RUNS 

Fee g 

٣٢۲۱7۶1۴۴۲۱۱۱۲۶۰ ۱1 50۴ ۲ ۶ 
۲0 

FL8=1 EVALUATES ٣٢" +7 77٦1 
FL8=1 

EVALUATE SPEED OF SOUND CONSTANTS 
Al-SQRTCGAMMALIXRIXTI) 

AG=SQRT CGAMMAG*¥RGXTS ) 

1 ,ب۵۸1 

۵۸41 0 - ٤۶ 


- о е ц ц ж тю «ж» тю ён «ше ён ён ён ён ён ён mE цю цю ца ца ц s s s s s s s ца ца цю «шт s s s s s s s s s s s s s s s s s s s s 
=< — ж — — — یب‎ ш= «ша «ае «ж» «ж» «а что «ж»  — — — аа" «а -«-! «а «а т өш» еш» ہت‎ г әш» шш «ж» тш» з یب‎ Чы w w q — ён w ہی‎ Фи s m s s s s s «ә 


EVALUATE CONSTANTS USED COMMONLY USED 
Cl1-CGAMMA1*12/C2XGAMHMA12) 
61026 6/272 
С5-(-2ХСАММА42/(САММА6-1) 
C6=(GAMNMAG-1)702*GAMMAG ) 
C7=(GAMMNAG+1 70 4*GANMAG ) 
C8=(GAMNAI-1)70C2*XGAMNAI DX ( 
С9-( САММА1 + 1) / (2Х( ОСАМНА1 - 1 
С10- ( САММА1 / (САММАТ - 1) ) 


Net e 


eee ee ee ён ч шш s s s s s s s s қыз қыз қыз қыз қыз қыз т” қыз қыз қыз қыз «шт өше қыза қыз қыз «шт ہے"‎ шт 
— — — — — — — — — — — < — — — = = = = = = – ہے — — — — – = = = = – – = = = ہے‎ шш шш шр — — — — — — =< — — — — — — — шшр == 


THIS LOOP INCREMENTS ۲۳۴۴55۴ ۲ ۶" ۶  >))")/+؟‎ ۱)۱ 
USES INTERVAL HALVING TO DETERMINE M3. 


рО 1 Ј=1,ММ, 1 

СОН5Т-РАР1(,2 

CALL INTHVCFCH,X1,X%2,XR,XTOL,FTOL,NLIM, 1) 
M3CJI=XR 
C3(J)=CC3 
۲۹9) ( (9 
X1=XR 
Х2=6.0 


С THE ROOT HAS BEEN FOUND AT XT ТР [QU P P ISO P. 


1 


Cat) GC CD OOO O OO fee, 


OOO 


K=J+1 
РАР| (К) -РАРТС .)) + А 
CONTINUE 


THIS LOOP USES P47P1 & M3 TO DETERMINE TZ TI ОЗЕ 


DO 5 1211 
FIND PZZPI CEQN S557) 


P2P1(1)=P6P1(1)ZC6(I) 
МСТ) =$50КТ(2х(Р2Р] (Т)+(1/Р2Р1(1))-2)/ ССАММА1 *  САМИАТ-1 2х 


&((GAHHA1+11)Z(GAHHA1-1)TP2P101) 70 


FIND MS & SPEED OF SHOCK SC Etta asi 

МУ 2-50КТ ССР 2Р А а Е ОСКИ 

US(I)=MS(I)%xA1 

FIND P37P4 

Р5Р4(1)+С5(1) ХХС5 

РаРУ(1)< | /РЗР4(Т) 

EVALUATE MAXIMUM DURATION OF UNIFORM FLOW IN THE REGION 
BETWEEN THE SHOCK WAVE AND THE CONTACT SURFACE. 
REFERENCE PAGE 8641 


T2(I)=(2Z(P3P&6(I)%XxC7)-(6xA6Gx(1-FP3Pa(I)xXxC6)x(1ZP3P6(I))%xC7)Z/ 
(GAMHAG-1)XAlXSQRTCItCGAMMA1412)XCP2P1CIO-1)07 
2XCGAMMA1))))0XL6/A6 


T3 IS THE THEORITICAL MAXIMUM DURATION OF UNIFORM FLOW 
ІН THE REGION BETHEEN THE CONTACT SURFACE AND THE TAIL OF 
CENTERED EXPANSION HAVE CREGION 3). (PG 433 PBLM 8-6) 


Т5(12:(14/А46)%((САММА6-12/22Х%((РАР3(1)%%(3-САММАЯ6)) 


&X*¥(1/(GANHA4X4)))/(1-(P3P4(1)XXC6)) 


TC & XC ARE THE TIME AND DISTANCE AT WHICH THE WAVE FORMED BY 
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FILE: SHDCK2 FORTRAH А! 


B ТИЕ КЕЕШЕСТ ИН ИШЕН НЕ CENTERED EXPANSION WAVE OVERTAKES THE 
C CONTACT SURFACE. CREF PG433 PBLM 8-7) 
С 

ТССТ) <( | 4/ Ад) Х2Х( РАРЗ(Т) ХХС7 ) 

ХССТ) < (| 44 / (САММАЯ9- 1)  (1-(РЗРа(Т) ХЖХС6 ) ZC PGP3C1)*x*C7))) 
C 
Дије - == -=_____ ____-_____- = - = == = = = ши ннишшенв=кмнњнншегесев-== = =—=-=-==<=-=-=—=========== 
E 
Ë CALCULATION OF REFLECTED REGION PARAMETERS ACCORDING TO REF 2 
С 
C VELOCITY OF CONTACT SURFACE УС 

УМССІ) =2хА1/ ССАММА1+1 )хСМ5(С1)-1/иМ5(І) ) 

Е PRESSURE AND TEMPERATURE IN REFLECTED REGION 5 CEQH 2.31 : 


Р5РІ СІ 2 а (ОеХБАННА 12845 С І ) жх2- ССАПНАЇ - 122/ ССАМНАЇ 51) ) 
&ACCCSXGAMMA 1-1 ХМУСТ)  Х2-2Х (С САНМА1-1))/ ( (САМИА1-1) “М5(СТ) Хх2+2) ) 
Т5Т1 (СГ) <( (2Ж(ССАНИА | - 1) ХМЪСГ) ХХ 4 ( 5-САМНА| )) 
&x( (ЗХОАПМНАТ- 1) ХМ5СТ) ХЖ2-2Х( САММАТ-1))/ ( (САИМАТ+ 1) ХХ2ХИЗ(Т) ХУ2)) 

С САШЕШЕАТЕ VELOCITY OF REFLECTED БНОСК ОК СЕОМ 2.58) 

КВО РОУ СТО КОСА РОВАСБАННА ۶ С 190 )/ 
&ССОАММАТ+ УУ ССАМНМА1-1)-17Р2Р1(С1Т))) 

e OBSERVATION TIME AND DISTANCE INREFLECTED SHOCK 
DELTINCID- LLISCUSC CIT) -YCIA CUSCINI СОКО НУ СК І))) 
DISTCID-DELTIHMCIOXURCI) 

8 TIME TO EXPANSION WAVE REACH CONTACT SURFACE CEQH 4.7) 

СЕ( ( САММАЯ+1 )/ (2Х(САММА4- | )) ) 
ТСССТ) + 2Хх| 4Х( 1 + (САММАЗ- 1) ХМ5(Т)/2) “ХС/АЯ 
Pe LCT eC C1 AC C2xXAEX OSC 10> 42) OZ (CGANMAI+I DO XNSCI)D)) 
Е ПЕ ЕШТЕ BETHEEN CONTACT SURFACE AHD SHOCK 
DENISE OTCCCIOXCI-S2ZXCHSCIOXX2Z-197 6COADPBAI-tTI)XHMSCI)Xx2)) 
ppc ZP5PICIOZP2PTCIO 
5 CONTIHUE 


С SECOND METHOD OF CALCULATIONS 

X1-.1 

Х2=6.0 

ХТ0120. 0001 

FTOL=0.00001 

1-0 

НЕТМ=50 

XR=.5 

DO 2 J=1,MM,1 

COHSTzPGP1CJ) 

CALLSTIHDBHBECECHLI,XI X2 ХЕ ТОС РТО „МЕМ, ТІ) 
XHSCJ)zXR 
XP2P1CJ)2C2XGAMMAIXXMSCUJ)XX2-CGAMMA1-10)/7CGAMMA1*1) 
APSPáCJOozIZCPGPICJOZXP2PICJO) 
АКФОКТСОО-С(бАППА ОХХИ ЕЈ АОСОАММА] 1) 7 ХИ5(Ј) 77 2+2) 

ПЕТОЈ ИА ХАЈДЕ ЕМ 5۵ ۳1۲۱۸۱۱۰۱1077 27:794) ۴ل )2۶5×115 ص6۸1111۸1-1(7))‎ 
ЕО ССССАМНА ЕЮ ад ха> 


C XP2P1CJ)zC2*XGAMMA1»xHSCJ)Xx2-COGAHMA1-120)/CGAHITA1*1) 
۸۲۰ ۳)٥ у= (C ((GAHHATrT1)ZCGAPH4A 1-1) y+2-1⁄/ZXP2P102))2Z 
&C1+CGAMMAI+1 9/7 CGAMHAI-1)X1/XP2P1C J)) 
C XP3P4(J)=1/(P4P1(J)/XP2P1(J)) 
С ОСКЕ СОСЕМА ОУ ИЗ ОЭ 92 АСССАННАТ- 1) #45 (Ј )хХХ 2+2 ) 
E ХР) = С СБАНМАТ МЫС Л ХХ (БАМПНАТ-1)/2)) Х( (ОДММА1- 1) /2ХМ5( |) ЖЖ 
С &2*1)/CCCOGAMMA1*1)/2) XX2XMS CJ) XX2) 
X1=XR 
Х256.0 
1-0 
С 45 HAS BEEN DETERMINED BY INTERVAL HALVING 
C ПАПА ООШ ЕТОГУТУИИГІЦУ/ХЕ,5.0,0.0001,0,.00001,0,50/ 


ÜRUSCJOSURCJOZUSCJ) 
P5PGCJ)-ZP5P1C J)/PGP1CJ) 
е ۲00 1 
۲90 ۰ ۰ XxPG 
РТСЈ)гРаЛРАаРТСЈ) 
ОКОВЕ ТА 
ИЕШЕ А ОНИ) РАРЕСЈО, РаСЈ)», ОЕСЈ), БЕСТТМСЈ ) 
г CONTINUE 


С 129 


FICE: 


C 


СОШ 
С 


602 
С201 
С 

509 


С? СУСНОЛОМОСОУ СО) 


SHOCK2 ҒОКТКАН А1 


зар із іе «е «т і сш, ہے‎ «і сг г г са еш, «ше «шә «ж. «ше ош, سی‎ еш. «ын іш” «ше іш «а т тш эш жж» іш” тш сы, ыа ажын шш сш, шш шт т=з «шт «же ош, сш, шш тш ше -4- шш” шш” «ше «ше ہہ سی‎ быт «ше шэ «ше “ше сш» 
се «з с «ы іш еш» «ше «ше «ше «шт чи очен سج‎ — — — «к. «а еш» «ж» «ж» «те «ж» «ж» т «шт «ж» «ж» «шт шт аш” «шт шт” یہ‎ «сц і” «ае еш» «ж» «ж» «т іт г «т іт с іт і” г з” з «ж» «ж» «ж» = = = = = = 


IF СРЕЗ. ЕО, ООО 
THIS SECTION DETERMINES F5 T15 AHD DETERMINES TESTIME ASTA 
FUNCTION OF NOZZLE THROAT DIAMETER. 
НН=1 
ро 5 1-1, ММ 
REFLECTED MO AVAILABLE TO DRIVE NOZZLE 
СЁЗВЕТ(САММАІ/ЕК127Х(1-(бАММА1-12/2)хх(-1хС9) 
THASS(I)=P5(I)%x(UR(I)xXDELTIM(1))x5.1617xDT1xx2ZC(R1xT5(1) oa 
VIRITECX,X) РЗРІС(І1),Р5013, ٔ 1 ۵( 
CALL ۱۱822۲٤۲/٣ 1( ۰۱5 )٦٢3 NS TMA e CID, CFLS) 
FORMAT('P6P1=!,D9.35,2X, 'P1=',D9.3,2X, 'TOTAL MASS 10۰ ,109 , ر2 ر9‎ 
&'КУН =',15/) 
[1 
CONTINUE 
DO 309 I=1,100 
WRITES, 602) DINRCT TCI ITT ETTI SIC) 
ҒОРМАЛОВ СЕ 5/20) 
FORMATC'P5=',EËE9.2,2X 15 Ее ала ۹ 
&, 2X, ‘TOTAL MASS ٔ)۶ ۵( 
CONTINUE 
НЕІТЕ(8, 605) 
FORMATC/,'THROAT DIAM',2X, 'RUH1',6X, 'RUH2',7X, 'RUHS',7X, ' RUIHG',7) 
HRITEC8,601) )۲۹۲ ۱) ٣7 ,۲1۱)17 ۰ ٦ 


с, аз іт іт ہت ہت‎ «т іт і” г” іг і” і” і” ц” іс” і” і” i e 
- — — — — سے ہہ — = — = = – = — — — سے‎ — — — — — — — — — — — — – = == = – – – — — — == — — жш سے‎ == жш — — — — — سے‎ — — — — шю کے‎ 


бш» іт «ш» =з «шә «шы‏ ہے ہت Чыт - - «ше «м «ше «ше «ше‏ سث ہت ہے іст‏ ہے ہہ s‏ ہہ іт т «ц — іт у іт «а s s «ц s s а s s іт‏ ہے г ст іт‏ - سے سے یہ 
=з =з ==‏ ہت «ше «шт‏ ہے ہت «ше еш» т т «шт «шт «шт‏ ہت «те‏ ہت بہت ہب «ше «ж» «ше еше і” «т еш» еш» жж «м. «ж» «ше «ше «ше «же «ж»‏ سب ہہ ےہ ہہب т «шт‏ سب ہے سب سب «шт‏ یہہ че аа «а «а «т «т «г «же «же «ж» = т‏ 


1.۴۰.۰۲) ۲۲ 7۶ء۰‎ 
THIS SECTION DETERMINES PS TS AND DETERMINES ٦۱ 6510 НЕ 
OF PT/PTO FOR A GIVEN THROAT DIAMETER 
N= 
DO 82 1-1,ММ 
REFLECTED MO AVAILABLE TO DRIVE NOZZLE 
С=5ОКТССАММАТАЕТ О = ССАИМАТЕТ 2С Ор 
ТМА55( 1)-Р5СІ1)ХСИУВСІ)ХрРЕЇ ТІМС:(1))х5 авг река аса ٦ 
CALL NOZLEZCPSC1), 1501); Nines Cee 
МНЕНН+ 1 
СОНТТНУЕ 
DO 89 І-1,ММ 
WMRITEC8,202) PTPTO(I),FT(I,1),FT(1 ет o ICT O) 
РОКШАТСЪКЕЗ: 5 ака 
CONTINUE 
HRITE(8, 803) (P4P1(I),1=1,MM) 
FORMAT(Z ,8C('Fa6apl=' ro 502770 


= = =з ہت ہت == ہہ‎ шш «ш» > Eccl 
© — سب‎ «шт = = ~ = — — — — — — — — — — — — — ہے — — سے‎ быз «ш» «ше «ше ےہ‎ сш ш ж «ше «ше ہے‎ = == == = === == = = = = = = = = = = 


IF ٠٠۴ ۲2۰6 6۹۶ 
DETERMINE TEST SECTION MACH ) ")ٴ۶‎ ۷۳ 


CALL EXPANCTHR1,DT1) 
SUBROUTINE EXPANCTHR1,DT1) 
REAL X1,X2,XTOL, FTOL,DT1, THR1 
INTEGER ٦ 
DIMEHSIOM MTC2000,DTESTC2OUUOSDDDIHRCSS 
СОМИОН УА/ CONST, GAMMA], GAMMA4, AGA1SQ,C1,C1S5Q,CC3,CC4,C5,AGA1 
COMMON 7C7 ATASTR,MT,DTEST,C9,DTDTHR 
HRITE(8,300) 
DIEST(1)=DT1 
X1=.10 
XR=50 
DO 30 J=1,8 
Х2=100 
XTOL=0.0001 
FTOL=0.00001 
120 
HLIM-50 
АТАЗТЕ<( ОТЕЗТС 1 372 
MRITECX,X) ТНК1, РТЕ5ЗТС 1), АТАЗТЕ, С9 
CALL INTHVCFCN3,X1,X2,XR,XTOL,FTOL,NLIM, 1) 
MTC JI=XR 
X1=XR 
X2=100 


K=J+] JE SNC 


FILE: 07. FORTRAN АІ 


ПЕСКЕ ОЕ Э ТИСИ 
0٢٣0٢٢۲١ *>(ل‎ 0۲31) 71 
БТ =; HELP" 
ДПЕ СО БОГ ОНА У ОТОН СОЈ У DTEST( J), ATASIR,THR1, DT] 
300 FORMATO AA TEST", 1X," DTEST/DTUBE", 1X," DIEST!, 2X, "ATASTR' 
а. 9х, К" рТОВЕ" ) 
301 PORMAICGCF9>3),202X%,E9.3)) 
30 CONTINUE 


№55 ТЕГЕСЕ E070)GD TD 105 
C IIS ECTION EVALUATES TEST CONDTIONS FOR VARIOUS 
C TEST CHAMBER AREAS GIVEN A FIXED NOZZLE THROAT DIAMETER 
C OUTER LOOP IS FOR PGP1 CHANGE, INNER LOOP EVALUATE CH AREA 
DO 22 J-1,MM 
HRITEC8,400) 
DO 21 I=1,8 
ИЕС x) 0 
C шары да ۸۰۱٣٢٣١٦ ۰۲) ((,۲ 0)۱ ۵٣ 
С406 ПАСА АНКЕ GX "РБ =t 6X; TFPI"; 5X "T5! 8X, TIT! 5X; 'TSTART') 


C405 ۲۱۲۹۸۲۷۱) (ررے.809)‎ 
САИС сие то TSC STSCI, PICT LiCl), ISTARTICID, THR1,C10) 
T CIOS CII POPI CID, PICI) 


І 
ӘБЕС НОЮ НС І),Р5С(0),РТС1),1500),1 
8, ТЪТАКТ(Т ) 
400 ОО ТЕГ ACR SSA, Р5', ах, 'ТЕЗТ Р', 9х 
сале сацы OS, MS rax, "POaPI',6X, 'Pl',7X, 15 
401 БӘКЕ СЕС 22 1Х),Е6-. 1) 
21 CONTINUE 


С DUMP CHAMBER SIZING EVALUATION 
МКІТЕ(8,411) 


е 
” 


о ах, 
ТАКТ") 


кіс Ба ра 


XR = 50 
X2=100000 
XTOL=.0001 
FTOL=.000001 
1-0 
HLIM-200 
С ORTTECS,; XD) ААС11),РОТСІ1) 
CALL INTHVCFCHG,X1,X2,XR,XTOL,FTOL,NLIM, 1) 
VDVOCI1)=XR 
C MBITECS,%) AACI1),PDICI1),VDVOC11) 
X1=XR 
X2=10 
VDCIIYSVDVOCIAOSCGURCJOCDELETTHCJOXS. I1817XDTI*x276) 
RITECG GIZLI FSC, TSC, UDVOCII1J;VDCI1),AAC11) 
ШЕН ов, Эх, Vor, ok,  VDVO'T,7X,°VD"', 8X, 'PTPIO', 
) 


&6X 
412 ҒОКМАТ( 5(Е9.2,1Х)) 
460 CONTINUE 
461 CONTINUE 
ВЕТРЕ, 9022) DELTIMNC J), 10C1,J) 
402 FORATE TEST TIME OF SHOCK TUBE DI. 


CIDE TOR SPS OUT OF NOZZLE =',D9.6) 
С So vie TO STARIAROZZLE CEMPERICAL ESLIMATE)',D9.6) 
гг CONTINUE 


зе ц «ы шш с «ы з іа аш" «а «ж» «ж» «ше «ж» «ше «ше қы» «ы т ст қы» қы» қы» қы» қы» қы» қы» «м» se ты» қы» қы» «ж» қы» қы» қы» жы; ты» «ы ц” сь سب ہے‎ қы» қы» | ew we SK шш з ы ہے‎ ты» уш س س‎ Re үш 
— سے سے ہے ہے‎ чш سے‎ «ме «ше «м» «ше. «ше ہے‎ «шт аж. «ме қы» ат «о іш -" = а іш іа і а іт «а і і тт іш іш == == = шш = سے ہے‎  — «ае ہے ہے‎ шш «ше ہے ہے = سے ہے ہہ ہہ ہے‎ «ца سے — ہے ہے ہے ہے‎ 


105 K=MM 
IF (FL6.EQ.0)GO TO 699 
ИВІТЕ(8,695) 
WRITE(8,696) (PGP1(JS),MSCJ),USCJ), VCCI), URUSCJ),T5CJ) 
STINK) 
5 FORMAT(2X, 'XP4P1',7X, 'MS!,9X, 'VS'! ,9X, ' VC', 9X, 'URUS', 
&5X, 'P5', 8X, 'T5!) 
696 FORMATC6CF9.5,2X), E9. 3) 
699 IF (FL1.EQ.0)GO TO 108 
НЕІТЕ(8, 201) Тер! 


FILE: ی۶‎ 2 FORTRAN Al 


Ceo OA Oi) GGG) 


OO 


С СУСС С 


OO 


PFIRITE(8,200) (PáP102),M5SCJ); XPoP10J);XP5PoO2CUJ) ۹۹ ۵٥ 
&,XRO2R1(J),XT2T1(J),XH5S(J),URUSC(J),J=1,K) 
200 FORHATCOC(F11.5 2 ٦ 
201 FORHATC2X;, 'XPaPI'S7X, ۶ی ×8 ر "45ا۶‎ ×۲۲) СХ ИО РЕ 
۱ھ ۴ ر ڈ4 ,۲۴0027۳0۷۱1 ر6 , 8ش‎ 2711٣ ОХ БРАК Ы 
НКІТЕСЕ, 205) (PPID SMSC МАЕ МН Goo 'ٔ ер 4۵0< 
&, P5PáCJO, DELTIMCJO, J-1, KD 


205 FORMAT(7(D9.6,2X)) 
108 НЕІТЕ(8, 204) 
МЕІТЕ(8,205) К1,Т1,САММА1,К4,Т4,САММАЯ4 
206 FORHATCSX, 'PGPI'8X, MS ЕХ ТЕШИ ОИУУ ЕВА 


&,7X, 'P5SP6',5x,  'D I TI jr > 
205 ҒОЕМАТСТКІ-",Ғ-9,2,2Х,7Т71-”"/ ЕО А ОВИ ЕЕ УСО 
а7,'К4=',29,2,2Х,'19=',29.2,2Х, 'САММА8=',29.2) 


ст «ны تہ ےت ےہ‎ «к» өн» «т» «ше өн» ен» «не «ме ыз «ше қы» «м» «ғы «мн г т ц. пе «ғы «МЫ өт» ہج‎ с ہت‎ тж гц цс жш сь цц өт» «Мн цю шт жм» жы» жы» 
<< че ңқ ғ» ж. ہے‎ - а. ы. ъ= оь «ы ж. «ы ہے‎ с» «т ас - $ «ж» ھ‎ сь се се ہے سط‎ әже «ж» «ше - $ $ $ ат, еш» сы» «ны «шт «шт «ше «ж» «ж» 


ТЕ (FL7.EQ.0) GOTO 109 


OUTPUT DATA IO A FILE 
WRITEC8,103) 
І1РІТЕС8,100) СРАР1СІ), МУСІ), М2(1), МЗСІ),АВА1, 
8۲5۶۲[1)12( 71511٤120 7٤د6‎ )]11( : ٣٣1 ٣٦٦3٣٠٣٢ ٴ6‎ ٦.7 
МКІТЕ(8, 106) 
НКІТЕСЕ, 197) «(РеРІСТ»,ЮКСІ), РІ ЗТ ОПОБеХ СОТ му ПССОГО ВЕШ ЕЕ 
&DELTA(I), 21), КО) 
106 FORHATCGA, 'PGP1",2xX,*VEL REEE 39%, DIS PEr D OS 
& "EXP POS',4X,*EXP TINE’, 4%," DE О 1٦١ 
ИКТТЕ 8, 1083 
107 ҒОКМАТ(8(Ғ9.4,2Х)0 
PIRITE(8,105) (P6P1CI),1T2C1);,T5 1 7+ ITE) 
100 FORHAT(11(F9.6,1X)) 
102 FORHMATCaX, 'PGPl',8X S TZ OX, TFS OS е 
105 FORHMHAT(6X,'P6P1'!,6X, "HS", 
SIX, ۲۲۲2٢ ۲۷۹ھ 6 ر‎ АЧА X ۹ ۳ 
SOX TSI Ox, 
а1Хх, !РЪРа!,бх, "Р2Р|!, 6х, ВоВ ОХЕ 0٦ 
104 РОБМАТС2Х, "РОРІ"',12Х,'Т218Х, "ТР 5 ЗОНУ 175 
105 РОКМАТС5(Е10.6,2Х)) 
109 Х=1 
STOR 
END 
REAL FUNCTION FCHCX) 
REAL X 
СОМНОН/А/ CONST, GAMMA], GAMMAG, AGA1SQ,C1,C1SQ,CC3,CC4,C5,AGA1,C9 
CC3=(1+(GAHIHIAG-1)XXZ2) 
CC&%&=(CC3xx(2xGAHHA6))xXx(1⁄/Z(GAHHA6-1)) 
FCH2CI*C. 5BXAGATSQOXCCGAMMA1XX/CC3)XX2) 
ЯХ(С1 +5ОЕТС ( С150+94/ (АдАТ54Х( САММА1 ХХ / ССЗ) ЖЖ2)))) ) 
&¥CC4-CONST 
НКІТЕС8, 200) РЕСМ, бАММАІ, GAMMAG, AGA1SQ,CC3,CC4,CONST 
200 РОКМАТ (7 (Е8 .2)) 
КЕТОКН 
END 
REAL FUNCTION ЕСН1СУ) 
REAL Y 


COMMON/A’ COHST, GAMMA], GAMMAG, AGA1SQ,C1,C1SQ,CC3,CC4,C5,A4GA1,C9 
ЕСНІ з (СОХОАННА 1 ху хх2-АМНА1 41 7 ( GAMMA1L+1))*C1-CCGAMMNAG-1)7CAGAI1 
&XCGAHHATT13)x(Y-1/7Y)))2XxC5-CONST 
۲۴۲۱1٢۶ 1ء ۴6۲۱۱۷۰6۴0۱۱57 ۰ر‎ CS 
RETURN 
END 


REAL FUHCTION FCN3(Z) 
REAL Z 


COMMON /С2 ATASTR,MT,DTEST 
COMMON /А/ 601151 , 6۸۲1۱4۸1, ,6۸۱4۷۸۹4 ۸4۸159 ۴1 ۴159,6۴) ۴)4 ء۲‎ ۹ 


32 


FILE: ЭШШЕКС2 ҒОКТКАН А1 


С HRITECX,X) ATASTR,GAMMA1,Z,C9 
FCHN3:ATASTR-1/ZXC2XCIT CGAHMA1-12/2xZxx2) 
%/ СССАМНА 1 11 )))ж2С9 
HWRITECX,Xx) ATASTR,GAMMA1,Z,FCHS 
RETURH 
END 


REAL FUHCTION FCN4CHD 
REAL H 
С 6ئ۲۲‎ ٘ 11 
DIMENSION P5€200),7T5€200),AAC8),PDIC8),MTC 200) 
COMMON/C^/ ATASTR,MT,DTEST 
СОННИ ВР, 15,ТЕ, АА, РОТ, ‚11, 112 
FCH&a=l-(1+PDI(11)xXxT5(J)%wH/(P50J)xXT1)-AA(11)xx.716)xT1% 
з ої нич ]00)x=Xx35.5SZ(AACI1)xT5(J)) 
НЕВЕ, ГИ AA(11),PD1(11),P50J),M,T50(J),J,I1 
вы ое, PDI-tyj;E9.2, 'P5-',E9.2, 
НЕ ТА ЕӨ 2: ІХ, Ғ6.1,1Х:Ғ4.1) 
RETURN 
END 


DIDIS 


SUBROUTINE NOZZLECX,Y,NN,MTOTAL,CONST2,FL3) 
THIS ROUTINE DETERMINES TINE FOR NOZZLE FLO}! FOR A GIVEN THROAT 
DIAMETER. 
REAL x8 X,Y,MTOTAL,CONST2,PI 
ГЕТЕБЕКЖІШУГІ 5 
DIMENSION ОТНЕС200),1(200,94) 
COMMON /А/ СОНМ5Т,САММА1,САММА4,А4А150,С1,С150,СС5,СС4,С5,АЯ4А1,С9 
С COMMON /B/ P5,T5,R1,HM,C8,DTHR,DTHRDT,T 

COMMON ^B/ R1,MM,C8,DTHR,T 

) 8  ّي6‎ 

ОТНКС1) <, 0065 

ро 0 71 

ИЕДИ ООУ ОКИ ОСТ  SXXCS I1 (CONSTZ2*XFIXDTHRCI)XX27/4X%X 

&x(GAHHA1-1)) 
J*I*1 
DTHRCJOZDTHRCIO*.0005 

С МЕТИ ОВ 200) DTHRCID; TCI, NN) 
C200 гокшашкесЕз. 5, ех)) 
C201 OFA ES Ев ах: БЕ", ЕЗ 2а иие тт 
C ое ИА ASS РОТ, Ед, г) 
50 CONTINUE 


OO OO 


É WRITEC8,201) X,Y,NN,MTOTAL 
RETURN 
EHD 
D 
C 
SUBROUTIHE HOZLE2CX,Y,HH,MTOTAL,CONST2,FL3) 
С ME SROINE DETERMINES TIME FOR NOZZLE FLOM FOR A GIVEN THROAT 
С DIAMETER WITH VARATION OF PT/PTO 


REAL x8 X,Y,MTOTAL,CONST2,PI 
INTEGER HH,FL3 
DIMEHSION DTHRC20020,T€200,4D0,PTPTO(C10D,FT( 6,10) 
COMMON 7AZ CONST, GAMMA], GAMHAG,AGAISQ,C1,C1SQ,CC3,CC49,C5,A4A1,C9 
C COMMOH ^B/ P5,T5,R1,HMHM,C8, DTHR, DTHRDT, T 
COMMON ^/B/ R1,HM,C8,DTHR,T 
601410۳" 7/۴7 6۲6 ۴ۃ‎ 
۲16 
рТНК(123.0065 
ФІЙЕІДО Ы 2 
DO 80 I=1,4 
PRC IHH) SHIOTALXDSQRT(Y)X2%(1ZPTPTOCI)x%xxC8-1) 
&Z(COIIST2xPIxDTHR(1)xx2Z6xX 
Х( САММАТ-1) ) 
J= I+] 
ВТС) =РТРТОСТ)+.2 
80 CONTINUE 


RETURN 
END 
С 
С SUBROUTINE EXPANCTHR1,DT1) 
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) ٤ 767 FORTRAN Al 


С›УСУСУСУСУСУСУСУСУСУСУСУСУСУСУСУСУСУСУСУСУСУСУСУСУСУСУСУСУСУСУСУ 


40 
40 


G G GG G СОИС Gy Gy Gy G) GG‏ ہہ یہ ہہ تی 


5 
3 


REAL X1,X2,XTOL,FTOL,DT1,THRI 
INTEGER I,WNLIM,J 
DIMENSION MT¢200),DTEST(200),DTDTHRC5) 
COMMON “AZ CONST,GAMMA], GANMAG,AGA1SQ,C1,C1SQ,CC3,CC4,C5,AGA1,C9 
COMMON 7C7 ATASTR,MT,DTEST,C9,DIDTHR 
06 0717 
X1=1.0 
XR=2 
ро 12 
XE 10 
XTOL=0.001 
FTOL=0.0001 
o 
HLIM-50 
ATASTREZCDTESTCJO/THRI12Xx2 
HRITECX,X) THRI,DTESTCJO,ATASTR 
CALL. IHTHVCFCHS,X1,X2,XR,XTOL,FTOL,HLIM,I) 
MT(I)=XR 
X1=XR 
K=J+1 
ВТЕЭПСКЛ- СТЕ ӘР 
ОТОТНКС .)) 751 
МЕТТЕС&,501 ٦٦) ۰۱ +٤ + + 70 
РОБМАТС2Х, 'МТЕЗТ* „2х, "ОТЕЛИ ВОВЕ С УШЕУ 
НКІТЕС(8, 500) 
FOPHATC5GF9 205 
CONTINUE 
RETURH 
END 


SUBROUTINE TCORD(Z1,F5,15,P1,11, ТОЕ CID, 
REAL х 8 BSS TS Ш ТИЛЕДИ ГНЕ 0 
COHMNON ХА/ С015Т, САММА|, САММА4, АдА1 59, С1, С159, СС5, CC6, C5, AGAT, CO 
соннон /Е/ РрТЕН(г00), Т(200),МТ, рТНЕРТ, ММ, С8 
COMNON / BZ КЕ, М, Со. ٣ 
WRITE(X, 403) 

НЕІТЕ(Х, 404) 2-۶" HI 7 0 

ЧЕТТЕСЕ XX) CID 

A5-SQRTCGAHMA1*RIXT5) 

۲٢۰۲۹7۸ ) 14+) 6۸۲1۱1۸ 1-1 ؛ء , 2 ے2(‎  ۵[ 

ТТЕТӘУ/СІҒССАМНАІ-1)/2 ۵٣٦ 

B IS THE AHGLE OUT OF THE THROAT 30 DEGREES 

В-ТАН(15Х5.1416/180) 

ТІМЕХО 


TIME-.76 XTHRI/ZC2XBXABO*XCODSQUORTDOZIODSOLSTIAZZI) 
&tl/^/l15*CCDSURTCZI)OOXXxS-]1)*1725 РТ БЕТ ٣ 
FORHATC2X, MACH" ,6X,°P5 =*,6X5'FI", 54, 15' 36k, ' 11', 5X, Поа ан 
НЕТТЕ(Сж, 404) Z1,P5,PT,T5, TI Cl 6٤ 


FORMATC7 (D9. 209) 
КЕТПЕН 
ЕНО 


пи m «ы жы» سے‎ cuum «е се са тив сы” а а” с” са” «шы» Ұлы» «ы» سے‎ = = = = = = = о «шы» = = = у = «а а <= НА سب سب سیب‎ қожа «шы» «шы «шы «жы» «шы» жина سیب‎ ўве «ае «аг а а = еш” а? 


THIS ROUTINE HAS DEVELOPED BY PROF. MILLER DEPT AERONAUTICS 


SUBROUTINE INTHV 
THIS SUBROUTINE FINDS THE ROOTS OF A FUNCTION , 
F(X) = 0, BY INTERVAL ۷٣ 


PARAMETERS АКЕ : 
ЕСІ “FUNCTION THAT CONPUTES VALUES TORT CX .NMUST ВЕСрЕСЕ ЕЕЕ 
EXTERNAL IN CALLING PROGRAN.IT KWAS ONE ARGUMENT, X. 
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ЕРЕ ٰ ی۶"‎ 2 ҒОКТКАМ А1 


OO.) Gat? DODD DD 6(7 6 0۱77ی(‎ 0000070070 


Со 0100 


о 


ОО 


ССС О) 


20 


Al; X2 ан а аце са X F(X) MUST CHANGE SIGH AT THESE PTS 
XR RETURNS THE ROUT TO TBE MAIN PROGRAN. 
TUL TIO TOLERANCE VALUES FOR X; F(X) TO TERMINATE ITERATIONS. 
ІШІМ ПРАТ TO NUMBER OF ITERATIONS. 
1 -A SIGNAL FOR HOH ROUTINE TERMINATED. 

۷ MEETS TOLERANCE FOR VALUES: 

I: MEETS TOLERANCE FORF(X) 

ас WET EXCEEDED. 

ТЕЕ ВЕ НОЕ OPPOSITE TH SIGH TO F(X2) 


۲۱۱۳۶۱۲۰۱۱ ۶۰۰۹۱۷۶۲ 0۱۱۲۱۱٢1565۱٢۱۴0, ۲۲۲ VALUE OF I INDICATES WHETHER TO 
RINT EAC VALUE OR HOT.I120 MEANS PRINT ТНЕМОТ.НЕ.О0 MEANS DON'T. 


»4 سب s“ ца” «шы» «шы» «шы» «шы» «шы»‏ سے سے س «на‏ سے سے سے سے سے سے سے سے سے کے سے سے سے سے سے ше «ы» өше өн»‏ سے سے سے سے سے سے سے سے سے سے سے سے سے تب ہے سے سے سے س س س سے سے سے سے سے سے «н? іца?‏ سے ست سے «ше ўны‏ ست سے 


REAL FCN, Xl, X2: XR, XTOL,FTOL 
INTEGER HDTH,L, J 
КЕАЁ Е1,Е2,ЕК, ХЕКК 


ен» سے سے سے سے سے سے سے سے سے سے کے سے سے کات سے کے سے = = =« سے‎ «шае са «шат «шы» «шы» «шш іна іш? سے‎ u іна іча іа іа سس سے ست سے‎ s іа іа” і і а” а с” с” с” сн x سے سے‎ оғыз «лы» олы» «лы» «лы» «лы» «шы» «жы» سے‎ — Ы 


сан ыа ре?) DIFFER IN SIGN. 


EI ٔ ی۹‎ 70( 

пали геши ла) 

IF то Ве. O. O THEI 
БЕЛІ 2201 
RETURN 

END IF 


~ о о سے‎ 


аа шашы а ыны аб POIHTIS CONVERGING TO THE ROOT 


DO ۶۹۔20‎ 

7۲٠۰۰٠٠۰0 

[P = FCHK XR) 

ХЕКК < АВ5(Х1-Х2)/2.0 

ТЕТЕ O) THEN 
PRINT 199, J,XRZER 

ENDIF 


a s‏ < ہے سن سے سے سے سے ہے حتہ ЧЫЫР‏ سے سے س دہ ہت سس سس سس سے س ФИШ?‏ ج Фина‏ ست ست سے سے سا ست سے س ست سے Фина‏ سے «шы» «шы» «шы»‏ سے ست سے سے سے سے سے سے سے سے ست سے سے سے کس Фина‏ سے کے کے шша шша Чита e шана Чити‏ حم 


CHECK ON STOPPING CRITERIA 
55۰۹65 ۰۰6 ۷ب ھ‎ 


۲۶٢۰٠٢٦٢٠۰۶02 J; XR FR 
RETURNI 
ЕНО ТЕ 


IE ۸۶۹۰٣۹٠٣٠۰۰۰ 1 
[= 


PRINT 203, J,XR,FR 
RETURH 
ЕНО IF 


>{ س سس سے ج سے سے سے سس سس سے سے سے سے سے کے سے ست سس سے سے سسس سس سے کے ست سسس سسا سے سسا سسا سسا Á X‏ باد سے ست سے ست کت سے کی سے ست سے سے سے سے سے س س س س س س س کس س س ше «ше «ше іе «ае «е іна іза‏ سے 


СОМРОТЕ НЕХТ Х 


0۷) 0.0) ТНЕН 
X1=XR 
[USER 
ELSE 
X2=XR 
F2=FR 
ERD IF 
CONTINUE 
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FILE: 76 ҒОКТКАМ А1 


6 АБА о مہۓ جو ہے ہے ہے سے ہے‎ аа کے ہے ہے ہے یج ےج ہے ہے ہے رے سض ہے ہے ہے‎ үү - 
С 
Е HHEN LOOP IS HORMALLY COM LETED LIN TS EXCEEDED 
С 
I 
PRIHT 200, NLIM,XR,FR 
КЕТУКН 
с 
С те тшо ашы еше жш m tus оз мө же ше a е е сы се е == = == == کے بت بت بت‎ ОТ ОРНОГО РИСУ ОО М 
С 
199 FORMATC' AT ITERATION’,123,3X,"' X=t',E12.5,6X,' F(X)F И 
200 РОКНАТ(/! ТОСЕКАШСЕ НОТ МЕТ АРТЕК ",І4," ITERATIONS RK ZE 
1 Е2.5, ۶۰۲۲) ۶۹٤+ 0 
г01 FORMATC7' FUNCTION HAS THE SAME SIGN AT X1&X2") 


202 FORMAT(C7* TOLERANCE MET IN ',1%,' ITERATIONS X -',E12.5, 
1 F(X) > TE SD 
205 ЕОКМАТ(/' Е TOLERANCE MET IN  ',I14,' ITERATIONS X - ',El12.5; 
1 ۴ ۴)3 
ЕНО 
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